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Abstract The level statistics in the two dimensional disordered electron 
systems in magnetic fields (unitary ensemble) or in the presence of strong 
spin-orbit scattering (symplectic ensemble) are investigated at the Anderson 
transition points. The level spacing distribution functions P(s)'s are found to 
be independent of the system size or of the type of the potential distribution, 
suggesting the universality. They behave as s^ in the small s region in the 
former case, while s'^ rise is seen in the latter. 
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The random matrix theory, first developed by Wigner in 1950 's in the field of 
nuclear physics, is now widely applied to many problems in solid state physics 
such as the energy spectrum of fine particles, conductance fiuctuations in 
mesoscopic systems and the persistent current in mesoscopic rings In 
metals where the electronic states are extended, the level repulsion is strong 
and the distribution function P{s) for the level spacing s is well approximated 
by the Wigner surmise. On the other hand, in the insulating regime no 
level repulsion is expected since the electrons are locahzed and -P(s) becomes 
Poissonian. Natural question that arises from the above facts is how the 
properties of the level statistics are changed when the system undergoes the 
localization-delocalization transition, i.e., Anderson transition. 

To answer this question, much numerical effort has been devoted to the 
three dimensional (3D) Anderson model. Shklovskii et al. are the first to 
notice that P{s) at the Anderson transition becomes size-independent [Q], 
and the critical exponent v for the localization length can be extracted from 
the behavior of the size dependence of P{s) near the Anderson transition. 
Evangelou et al. P as well as Hofstetter and his coworkers made de- 
tailed analyses on the level statistics at the critical point, and P{s) is again 
shown to be size independent and to be different from the Wigner surmise 
or Poissonian. 

The analytical calculation of P{s) by Kravtsov et al. has suggested 
the possibility to obtain the critical exponent v from the asymptotic behavior 
of P{s) at the Anderson transition. They predicted that P{s) falls off as 

P(s) ~exp(-5s2-^), (1) 

with 7=1 — l/(iz/, d being the dimensionality of the system. It is further 
suggested that the whole behavior of P{s) might well be described by 

P(s) = v4s^exp(-5s2-^), (2) 

with A and B determined by the normalization condition, and /3 by the 
symmetry of the system. Setting z/ 1.4 we can estimate 7 at the three 
dimensional Anderson transition to be 0.76, and the exponent of the large 
s behavior of InP(s) becomes 1.24. The overall behavior of P{s) has been 
claimed to be well fitted to [^, Q 

AseM-Bs^-^^). (3) 
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It is further claimed that the same P{s) has been obtained in 3D disordered 
system without time reversal symmetry and that it is consistent with the 
recent study on the critical exponent in magnetic fields 

The linear rise of P{s) even without time reversal symmetry is contro- 
versial, since it is believed that the behavior of P{s) in the small s region 
is determined by the symmetry of the system, and in the absence of time 
reversal symmetry, P{s) should rise as ~ s^. The large s behavior of P{s) is 
recently questioned by Zharekeshev and Kramer 0, who have performed the 
large size diagonalization and asserted Poissonian behavior P{s) oc exp(— /ts) 
in large s limit, with k < 1. 

Thus the existence of the critical level statistics at the Anderson transition 
is well established, but it is fair to say that the functional form of P{s) is 
still controversial. 

In this confusing situation, it is very important to study the critical level 
statistics in other universality classes in the different dimensionality where 
the extensive studies on the Anderson transition have been performed. One 
of the examples to show the Anderson transition is the quantum Hall (QH) 
regime, where the Anderson transition takes place at the Landau band center 
with the divergence of the localization length as E being the energy 

measured from the Landau band center and i' ~ 2.4 [|10|, |ll|]. Independently 
of the above mentioned work on 3D transition, the level statistics in the 



quantum Hall regime has been studied by the present authors |T2|, |13[. It 
has been shown that the level statistics obeys the scaling behavior, and the 
size independent level distribution function P{s) has been obtained at the 
Landau band center. Another example is the two dimensional (2D) system 
with spin-orbit interactions which derive the localized states to the extended 
ones due to the anti- localization effect. 

In this paper, we discuss the critical level statistics in two dimensional 
systems, especially P{s), at the Anderson transition in QH regime (unitary 
case) as well as in the presence of strong spin-orbit scattering (symplectic 
case). We have obtained the energy spectrum by diagonalizing the Hamilto- 
nian with randomness for different systems and different sizes. The energy 
eigenvalues are then unfolded in order to make the average spacing unity 



|13| . From the whole spectrum we take out levels in a central region with 
a width of one tenth of the total width, and discuss their statistics. The 
number of samples for each system size is chosen so that the total number of 
levels for which the statistics are considered should not be less than 25,000 
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in the case of QH regime and 50,000 in the symplectic case, respectively. 

First we discuss the QH case. The disordered potential is assumed to be 
(5-correlated. Here we have adopted the random matrix Hamitonian which 



describes the energy spectrum in the QH regime [|TT], 0. The area of the 
sample is varied from 200 to 800 in units of 2^T^'^ with ^ the magnetic length. 
In Fig. 1, the P{s) at the Laudau band center is shown in both linear 
(Fig. 1(a)) and logarithmic scales (Fig. 1(b)). Different marks represent 
different system sizes. These figures indicate that the critical distribution is 



size independent, which is consistent with the scaling argument 



Now we apply (2) using the value u ~ 2.4 which results in 7 = 0.79. The 
result is shown by the broken lines. We see considerable deviation from the 
numerical data. Regarding 7 as a fitting parameter, very good agreement is 
obtained with 7 = 0.35, but this , together with 7 = 1 — l/dv, gives absurd 
value of z/(= 0.77) that violates the inequality u > 2/d |T1|. It should be 
noted that the small s behavior of P{s) is proportional to s^, consistent with 
the symmetry argument. 

A similar difficulty also appears in the analysis of P{s) at the Anderson 
transition in the presence of spin-orbit interaction. Using the Ando model 
[15] and estimating the critical disorder by the finite size scaling method 



15, 16], we have obtained the critical level statistics. In the numerical 



simulation, the ratio between the spin flip transfer V2 and the total transfer 



V = \JVi + V2 (Vi the non-spin-fiip transfer) is set to 0.5, and the critical 
disorder is found to be 5.75V^ and Q.TV for box distribution and Gaus- 
sian distribution of diagonal disorder, respectively. In the case of the box 
distribution the parameter W describing the strength of disorder is defined 
as the width of the box, while in the Gaussian distribution it is defined as 
2v^ times the r.m.s. of the fluctuation of the site energy. The factor 2a/3 
is introduced so that both distributions may have the same variance for the 
same value of W . The linear system size is 8, 12, 16 and 20 in units of the 
lattice constant. 

The results can be plotted in a similar way as in Fig. 1. In Fig. 2, P{s) 
at the Anderson transition in the presence of strong spin-orbit scatterings 
is plotted in the logarithmic scale. It shows that the critical distribution 
is independent of the system size and the type of disorder. The behavior 
of the P{s) is similar to the results obtained recently by Schweitzer and 
Zharekeshev [O] who have diagonalized very large systems in the case of box 



4 



type potential distribution. The critical distribution function P{s) cannot 
be fitted to the theoretically suggested form P{s) cx exp(— with 7 
= 0.82 (broken line) which corresponds to the critical exponent u = 2.8 
The best fit in the linear scale is obtained by putting 7 = 0.35. Though this 
value of 7 does not reproduce the critical exponent u, good agreement in the 
small s region supports s^-rise of P{s). 

In summary, we have obtained the level spacing distribution function in 
2D quantum Hall regime as well as in the presence of the spin-orbit inter- 
action. The level spacing functions P(s)'s at the Anderson transition are 
shown to be independent of the system size or of the type of the potential 
distributions. In the QH case, P{s) at the Anderson transition rises as 
while it rises as in the symplectic case. 

The numerically obtained data can not be fitted to the theoretically pro- 
posed functional form (2) with 7 = 1 — l/du. This should not be taken 
seriously, since the argument of ref. is valid only in the large s region, and 
the overall functional form, though very successful in 3D case, is only a con- 
jecture. The numerical data in the large s region qualitatively support the 
asymptotic behavior predicted by them. From the practical point of view, 
however, it is very difficult to determine quantatively the behavior of P{s) 
in this region, since the number of data corresponding to large spacing is 
extremely small. 

In 3D, P{s) at the critical points has been reported to rise linearly in 
s irrespectively of the symmetry of the system 0. In ref. 0, the time 
reversal symmetry is broken by gauge fields which can be absorbed into 
boundary conditions, and the system involves no real magnetic field. A 
study of more realistic systems will be necessary to clarify the behavior of 
the critical statistics in 3D unitary ensemble. 

The authors would like to thank Professor V.E. Kravtsov and Dr. T. 
Kawarabayashi for fruitful discussions. 
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Figure Caption 

Fig. 1 P{s) in the quantum Hall regime in the linear (a) and logarithmic (b) 
scales. The circles (o), squares (□), diamonds (o), and triangles (A) 
correspond to the system area 200, 400, 600, and 800 in units of 27r£^. 
It can be very well fitted to the formula As"^ exp(— indicated by 
the solid line. Applying (2) with 7 = 0.79, we end up with poor fitting 
(broken line). 

Fig. 2 P{s) in the presence of strong spin-orbit scatterings in the logarith- 
mic scale. The circles (o), squares (□), diamonds (<>), and triangles (A) 
correspond to the system size 8, 12, 16 and 20 in units of the lattice 
constant with the box type potential distribution, while the crosses 
(x) and pluses to the system size 16 and 20 with the Gaussian 
type potential distribution. The application of (2) with u — 2.8 and 
7 = 1 — is indicated by the broken line, which shows consider- 
able deviation in the small s region. The best fit in the linear scale is 
obtained by setting 7 — 0.35 (solid line). 
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